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^ ; Abstract 

t^j- ■ We demonstrate how the solution to an exterior Dirichlet boundary value 

problem of the axisymmetric, stationary Einstein equations can be found 
in terms of generalized solutions of the Backhand type. The proof that this 
generalization procedure is valid is given, which also proves conjectures 
about earlier representations of the gravitational field corresponding to 
rotating disks of dust in terms of Backhand type solutions. As a further 
result, we find that in contrast to the Laplace equation, arbitrary boundary 
values may not be prescribed. 



1 Introduction 

X 

Extraordinarily massive, compact astrophysical objects like neutron stars require 
a fully relativistic treatment. This motivates the study of the relativistic gravi- 
tational field of axisymmetric, stationary rotating bodies. 

The investigation of the exterior field, i.e. of arbitrary exterior boundary value 
problems of the axisymmetric stationary Einstein equations, can only be a first 
step in this direction, for neither the shape of the body (i.e. the boundary) nor 
particular boundary values may be prescribed. For a given equation of state 
characterizing the matter of the rotating body, the boundary and the values 
there result from transition conditions since the metric along with its first nor- 
mal derivative behave continuously at the boundary (in appropriate coordinates). 
Therefore, a subsequent step must be the combination of a procedure to compute 
the interior solution with the method to treat the exterior field in order to realize 
the transition conditions. 

Our aim is to show a way to handle the exterior field. In particular, we want 
to make use of the variety of exact analytic solutions that are available in the 
exterior region. Thanks to the facts that 
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1. the exterior vacuum field equations can be summarized in a single nonlinear 
complex differential equation - the so-called Ernst equation |], ||, 



2. the Ernst equation has been found to be the integrability condition of a 
corresponding linear matrix problem, 

it is possible to apply various soliton methods (see H3|-Hl3|) to create explicit 
exterior solutions, among them the Backlund transformations and the Riemann- 
Hilbert techniques. By extensive investigation of the latter, Neugebauer and 
Meinel succeeded in solving the boundary value problem that corresponds to an 
infinitesimal thin rigidly rotating disk of dust fl4], |15| , jlBf . In particular, they 
found their solution to belong to a 'hyperelliptic class' of solutionsfl7[], see also 



T8j-|21|. Further investigations as to whether this class of solutions can be used 
to solve more general boundary value problems corresponding to differentially 
rotating disks of dust have been carried out by Ansorg and Meinel[]22| and by 
Ansorg|23|. In the latter paper it was demonstrated that the Backlund type so- 
lutions [] suffice to approximate the exterior gravitational field of an arbitrary 
differentially rotating disk of dust. 

In this paper we investigate the question, as to whether the treatment of the 
boundary value problem of an arbitrary differentially rotating disk of dust by 
means of Backlund type solutions can also be applied to extended boundaries 
B, which can potentially form the shape of a rotating body. We denote by B a 
smooth spatial curve, at which regular boundary values are prescribed. The solu- 
bility of such Dirichlet boundary value problems for sufficiently weak relativistic 
boundary values has been proven by Reula|24||, see also J25], pBfl . In illustrative 



examples we will take spherical boundaries in order to make the formulas more 
transparent and therefore get better insight into the underlying mathematical 
structure. However, the general statements are valid for arbitrary extended B. 
The paper is organized as follows. At first, the metric tensor and the Ernst equa- 
tion are introduced. Then, solution techniques to solve boundary value problems 
of the real axisymmetric three-dimensional Laplace equation with the boundary 
B are discussed in order to prepare the relativistic treatment. These solutions 
will be given in terms of an analytic function H, which is defined on a curve 
T in the complex plane and satisfies H(X) = H(X) (with X 6 C and the bar 
denoting complex conjugation). Here, the curve T is closely connected with the 
boundary B. In this formulation, the solutions describe both a regular interior 
and exterior field, which (in general) assumes different values on the inner and 
outer side of B. The prescription of these interior and exterior boundary values 
uniquely determines the function H, i.e. this formulation permits the simultane- 
ous solution of an interior and exterior boundary value problem. 



lr The Backlund type solutions can be obtained by applying Backlund transformations to 
purely real 'seed'-solutions of the Weyl class. These solutions form a subclass of the hyperelliptic 
class, see[p3|. 
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The second section treats the hyperelliptic solutions of the Ernst equation and 
their generalization by means of a suitable limiting process. Similar to the formu- 
lation of solutions of the real Laplace equation, we find these generalized solutions 

1. to depend on an analytic function 7 defined on T (and, in general, not 
subject to further requirements) 

2. to permit the simultaneous solution of an interior and exterior Dirichlet 
boundary value problem of the Ernst equation. 

In the case of a flat interior field we find an explicit relation of the function 7 
and the values of the Ernst potential at the exterior part of the symmetry axis. 
In the section 3.1, we investigate the solutions of the Backhand type. Apart from 
regularity properties we discuss the generalization of these solutions which can be 
performed since they form a subclass of the hyperelliptic solutions. As a result we 



prove conjectures formulated in |[23|| , in particular that the Neugebauer-Meinel- 
solution can be written as a well defined limit of Backhand type solutions. 
Because of their mathematical simplicity (compared with the much more compli- 
cated hyperelliptic solutions), the Backhand type solutions play the fundamental 
role in our treatment for approximating solutions of an exterior boundary value 
problem by analytic solutions. This is carried out in section 3.2. As a result, we 
find that in contrast to the Laplace equation, arbitrary boundary values of the 
Dirichlet type may not be prescribed. 

In what follows, units are used in which the velocity of light as well as Newton's 
constant of gravitation are equal to 1. 

1.1 Metric tensor, Ernst equation, and boundary condi- 
tions 

The metric tensor for axisymmetric stationary and asymptotically flat vacuum 
space-times reads as follows in Weyl-Papapetrou-coordinates (p, (, <p, t): 

ds 2 = e" 2 VW + d( 2 ) + p 2 d^ 2 } - e 2U (dt + a dip) 2 . 

The field equations are equivalent to a single complex equation - the so-called 
Ernst equation 

(B/) A/ = (V/) 2 , (1) 

A-— -— — V-(— — 

dp 2 pdp <9£ 2 ' {dp' d( 

where the Ernst potential / is given by 

f = e 2U + ib with bc = a „ , b p = ar- (2) 

p p 
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The remaining function k can be calculated from the Ernst potential / by a line 
integral: 

k -f = (u P T-(urf + \e-± u [(b, P T-(b c T], 

k ~j = 2U p U c+ l -e-*% p b x 
In this paper we treat the following boundary value problem of the Ernst equation 

1. The boundary B is given by a smooth spatial curve which can be described 
by some positive analytic function r B : [—1, 1] — > M + : 

B = {(p, () : p = r^cost?) sin??, ( = r B (cos , d) cost?, $ G [0, 7r]} 

2. Along the boundary B we require regular boundary conditions of the Dirich- 
let type, i.e. / is given for (p, () G B. 

3. Regularity at the rotation axis p = is guaranteed by 

f(o,C) = o. 

4. At infinity, asymptotic flatness is realized by / — > 1. 

The treatment in section 2 will provide us a formulation / = /(j) which is suffi- 
ciently general to satisfy interior and exterior boundary conditions of the Ernst 
equation (|l|), i.e. for an appropriate prescription^ of (not necessarily coinciding) 
interior and exterior boundary values, there is a uniquely determined function 
7 : T -»• C [with r = {X G C : (\$X\ , dtX) G B}} such that the Ernst potential 
/ = /(t) is both regular within and without B and assumes the above bound- 
ary values at the inner and outer side of B. Since for physical reasons we are 
only interested in an exterior solution, the freedom in the choice of the interior 
boundary values will be used to restrict the function 7 such that the solution 
f — f(l) can De represented by Backlund type solutions. This is outlined in the 
third section. Furthermore, for the exterior solution we will impose the additional 
physical requirement of reflectional symmetry with respect to the plane £ = 
(see U?7[] ) which leads to 

r B {r) = r B (-r) for all r G [-1, 1] , 
f(p, -C) = f(j>,0 for all (p, C) outside 5. 
2 As already mentioned above, the Dirichlet boundary values cannot be chosen arbitrarily. 
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1.2 Boundary value problems of the Laplace equation 

In this section we prepare the general formulation / = 7(7) by considering the 
corresponding non-relativistic gravitational boundary value problem - i.e. the 
general boundary value problem of the axisymmetric three-dimensional Laplace 
equation 

AU = 0. 

Any real solution U which is is both regular inside and outside the boundary B 
can be written in the form 



with 



T = {X EC: (|9fX| ,3fcX) E B}, 

H : T — > C is an analytic function with H(X) = H(X) , 

W z = v/(X-C) 2 + p 2 , (4) 

< for (p, C) outside B or dt{X) < C 

> for (p, C) inside B and M(X) > ( ' 



The function H is uniquely determined by the interior and exterior boundary 
values of the potential U . For a spherical boundary (with radius re = 1) we may 
write H in the form 

00 

H(X) = ^(H^Xi- 1 + H^X-i) , H<f e R. 
3=1 

Then the coefficients can be read from an expansion of the boundary values 
in Legendre polynomials 

PJ X ) = -L^—Ux 2 - l) n ] 

of the cosine of the angular coordinate 1? with p = r sini?, ( = r cost?: 

00 

\imU(p,(;H) = J2 H j +)p j-i( C0S V) 

00 

limC/( P ,C;//) = -J2 H{ 3~ )p i-i( C0S #) 
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Discussion. 



1. Reflectional symmetry 

TT(n -t- U\ - I kU (p>('' H ) for (P,C) within , , 

lP ' *>' } \-kU( P ,(;H) for (p,C) without B ' l J 

with k 2 = 1 is obtained if tb{t) = tb(— r) and H(—X) = kH(X), i.e. an 
odd/even function H produces an even/odd exterior potential U in the 
coordinate (. 

2. The freedom of the choice of the interior field can be used to restrict the 
function H in an appropriate manner. For example we may choose dt(H) = 
const, for X e T, where this constant can be arbitrarily prescribed. For the 
above spherical case, this results in = ■ Then the function H is 
uniquely determined only by the exterior boundary values. Analogously, in 
the relativistic treatment (see section 3.2) we may choose the corresponding 
function 7 such that 3 (7) = const. 7^ for X G T, which ensures that the 
corresponding Ernst potential / can be approximated very well by solutions 
of the Backlund type. 



The hyperelliptic class of solutions and its 
generalization by a suitable limiting process 



were 



Meinel and Neugebauer [17]] as well as Korotkin |18|, |1J, see also |20], |2T 
able to construct the hyperelliptic class of solutions containing a finite num- 
ber of complex parameters and one arbitrary real solution to the axisymmetric 
three-dimensional Laplace equation. By investigation of corresponding Riemann- 
Hilbert-problems one finds, that this class can be generalized in an appropriate 
manner. The generalized class turns out to contain all Ernst potentials which 
are both regular within and without the boundary B, which assume (in general) 
different boundary values at the inner and outer side of B, and are sufficiently 
weakly relativistic, i.e. close to the flat space solution /' = 1. In this section we 
discuss these generalized solutions which allow us to approximate the solution to 
a simultaneous interior and exterior boundary value problem. 

2.1 The hyperelliptic class of solutions 



We adapt the hyperelliptic class of solutions as given in [TjJ to our purposes by 
writing: 

For a given integer p > 1, a set {Xi, . . . , X p } = {X u } p ^ of complex parameters, 
3 In the following, the notation {Xi, . . . ,X p } will be abbreviated by {X u } p . 
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and an analytic function H : T — > C, H(X) = H(X), the following expression 

X p dX 1 fH(X)X p dX 



( p r x{v) xprfx 1 r 



(6) 



with 



W 



\ 



(X + iz)(X - iz) f[(X -X V )(X-X„), z = p + iC, 



x ^ , rM dX 1 / #(X)X^' dX n ^ . 



satisfies the Ernst equation. The (z- dependent) values for the X^ v ' as well as the 
integration paths on a two-sheeted Riemann surface have to be taken from the 
solution to the Jacobian inversion problem (|7|). 

The Ernst potential / = f({X u } p ,H) [| is regular at (p, () = (\^s[X v ]\,^t[X u ]) as 
can be deduced from an appropriate combination of the equations (0) and corre- 
sponding rearrangements of the terms occurring in (|6|). In general, f({X u } p , H) 
assumes different values at the inner and outer sides of the boundary B, i.e. it 
possesses a jump along B. 



2.2 The generalization of the hyperelliptic class of solu- 
tions by a suitable limiting process 

This section contains a theorem which ensures that the above hyperelliptic solu- 
tions can be generalized in an appropriate manner such that simultaneous interior 
and exterior boundary value problems become soluble by means of these solutions. 
The theorem consists of three parts, the proofs of which are outlined. Moreover, 
we illustrate certain properties of the generalized class. 

Theorem. 

Given the analytic function 7 : V — > C which can be extended to some compact 
neighborhood G 7 about T, then 



1. for sufficiently small e, the Ernst potential 

f(ej)= \imf({X^} pi eH p ) 



For simplicity, we suppress the coordinate dependence and write only / = f({X u } p , H). 
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exists and is independent of the particular choice of the sequences {{Xl P 
and {H p }™ which serve to represent 7 by 

v 

7 (X) = lim H P (X) TT(X - Xl p) ) for X e G 7 . 
p— +00 -I- 

u=l 



2. the Ernst potential / = f{e r y) is both regular inside and outside the bound- 
ary B and assumes (in general) different boundary values at the inner and 
outer side of B. The exterior and interior Ernst potentials can be extended 
beyond the boundary B to the region B 7 = {(p, Q : £ ± ip € G 7 }. 

3. any sufficiently weak relativistic Ernst potential / which is both regular 
within and without the boundary B with (in general) different boundary 
values at the inner and outer side of B, can be written in the above man- 
ner as f — 7(7). The function 7 is uniquely determined by the Dirichlet 
boundary values of / at B. 



Proof. 

1. In order to prove the first one of the above statements we establish a relation 
between the function 7 and holomorphic functions a and (3 defined on G 7 
from which the Ernst potential / can be calculated via the solution of a 
Riemann-Hilbert-Problem. 



(a) To introduce a and /3, we follow the treatments of 14 , 28 and [25 
The Ernst equation is the integrability condition of the Linear Problem 



If + f)- 1 
if + f]- 1 



h a/; 

h A" 1 /" 
A" 1 /,, /,* 



where $ = <&(X,p,Q is a 2 x 2-matrix and the spectral parameter A 
is given by 



X(X) 



X-iz 



X + iz 

One explicit way to create solutions to the Ernst equation that are 
regular within and without B and possess a jump at B is to use the 
Riemann-Hilbert-techniques, by which we force the corresponding ma- 
trix $ to possess a multiplicative jump at T: 
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and moreover to be regular within and without T. The indices '+' 
and '— ' refer to the interior and exterior sides of T, respectively. The 
jump-matrix C = C(X) is independent of the coordinates (z,z) and 
can be cast into the form 

c(*) = (j$o i)' ^ = a(X) ' Wl 

If the functions a and (3 are prescribed, the Ernst potential at an arbi- 
trary point (p, C) can be determined by the solution of a linear integral 
equation, see [f28|j . 

As outlined in |29|] , the linear system can be integrated along the rota- 



tion axis, which yields in agreement with the above Rieman-Hilbert- 
Problem [p = 0, ( > -r B (-l)]: 



/ 1 \ / F(X) 

/ -1 J I 5(X) 1 

7 V ' (9) 

1 / \ / 1 B(X) N 

1 -/l F(X) 



The notation $ ± refers to the two sheets on which $ is defined; '±' 
means the sheet in which A = ±1 for p = 0. The functions F = F(X) 
and i? = B(X) f\ are defined in the complex plane and are regular 
within and without the contour r (F — > 1 and S — > as |X| — > oo) 
but possess a jump at T. This jump can be taken from the above 
Riemann-Hilbert-Problem: 

F + = aF_ , B + = aB_ + (3. (10) 

(b) We establish a connection between the relativistic parameter e as well 
as 7 and the functions F± and B±, now additionally depending on e. 
Thus, by virtue of (|10|), we get a = a(X; e; 7) and (3 = (3(X; e; 7). 

Given the function 7 which is holomorphic on G 7 , we denote by i3 7 
the Banach space of all holomorphic functions on G 7 with the norm 

||a||= sup for a & B 7 . 



We now define functions 



F± : [0, e 7 ] -> S 7 , S± : [0, e 7 ] -> B 



1 



D The function i? = B(X) is not to be confused with the boundary B. 
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that satisfy the following system of differential equations: 

= F ± (L± Sl ) + B ± (L±s 2 ) 



A 

de J 



e^F ± 



B l ± -\ 



(11) 



f^(L±s 2 ) + F ± (L ± s 3 ). 



Here, the Sj as well as the L± are introduced as follows: 
i. The Sj, Sj : [0, e 7 ] — > £> 7 , follow from F±,B± by 



l(F + 



F_)si + (B + + B_)s 2 ] tanh f = - F_ 



B-) 8l 



B\-l Bl-1 



S2 



tanh 



G 



G(B_ - fi 4 



s 3 = F+ 2 [2F + B +Sl + (£? 2 - 1) S2 ] 



G 2 



ii. Any function a e *B 7 can uniquely be written as the sum 



a 



with o" + and cr_ only possessing singularities without and within 
the curve T, respectively, and cr_ — ► as |X| — > 00. The linear 



operators L±,L± : £> 7 



Explicitly: 



1 

2vri 



X-F 

and these functions can be extended to 
By means of the substitutions 

1 + eF* 



£> 7 , extract these functions a±: 
L±a = =f<t ± . 



(L+a)(X) = -o+pO for x within r 
(L_a)(X) = a_(X) for X without T ' 



J 7 since a G £> 7 . 



F± 
B± 

S\ 
S3 



= 

the differential system (|TTT) reads 



e{B% - B*_)+e 2 s* 2 
Bl)+e 2 st 



d 
de 



( F*\ 

p* 

B* + 
\B- J 



( h!(e;F*,Fl,B* + ,B*_) \ 
h 2 (e;F*,Fl,B* + ,B*_) 
h 3 (e;F;,F*_,B* + ,B*_) 
V h±(e;F*,Fl,B* + ,B*_) ) 
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with the functions 

hj : [0, £ 7 ] x — > i3 7 

satisfying a Lipschitz condition with respect to all arguments in a 
sufficiently small interval [0, e 7 ]; the upper limit e 7 is defined in this 
manner. Note that for the continuity of the linear operators L± it 
is necessary to define them on £> 7 and not on the Banach space of 
functions analytic at T. 

As a consequence of the theorem by Picard and Lindeloff, the above 
system of differential equations has a solution which depends uniquely 
and continuously on the given initial conditions. This is the point at 
which we bring in the function 7: 

Fl(e = 0) = L ± ( 7 + 7*) , B* ± {e = 0) = -L±( 7 - 7*) , 
where the function 7* e £> 7 results from 7 by 

7*(X) =7P0 forXGG 7 . 

These initial conditions together with the differential equations (|Tl"|) 
yield F± and B± which only possess singularities without/within Y. 
Thus, the functions F and B obey a Riemann-Hilbert-Problem (|T0D , 
and the associated functions a = a(X;e;7) and (3 = f3(X;e;j) are 
uniquely determined. Consequently, we find a formulation / = f(e; 7) 
by solving the linear integral equation that has been mentioned in (a) 
for a = a(X; e; 7) and (3 = f3(X; e; 7).^ 

(c) If we prescribe 7 P in the form 

lp (X) = H P (X) f[(X - X«) 
i/=i 

then it can be shown that the Ernst potential following from the above 
reads 

f(e; lp ) = f({X^} p ,eH p ). 

The proof of this uses many of the solution methods that were devel- 
oped by Neugebauer and Meinel when they solved the boundary value 
problem of the rigidly rotating disk of dust. These methods are par- 
tially given in |p9| ; in full they will be treated in a subsequent paper. 
For any series {7 P }^ of the above functions 7 P with 7 P — > 7 as p — > 00 
and 7 P G £> 7 , the corresponding Ernst potentials / = f(e; j p ) con- 
verge for sufficiently small e since the functional / = f(e; 7) depends 
continuously on 7. Because of the above equality, this implies the 
convergence of f({X^} p , eH p ) as p — > 00. 

6 The solubility of this integral equation is ensured for sufficiently small e, see point 2. of 
the proof. 
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Remark: The differential system ([11]) can be explicitly solved if 7 is regular 
without T with 7 — > as | X\ — > 00. Then one obtains 

x^n 1 t 1 ^ sinhG o / ^sinhG 2 
a = coshG+ (7 + 7 )—2Q-> P = ~{l ~ 1 ) 2Q , G = 77 , 

from which the Ernst potential along the rotation axis can be taken directly, 
see ( fL3|) below. 

2. The proof of the second statement uses the relations a = a(X; e; 7) and (3 = 
j3{X\ e; 7). In particular, we have a, /3 G £> 7 for £ G [0, e 7 ] and moreover a = 
1 + 0(e) and /3 = 0(e). This means that for sufficiently small e, the linear 
integral equation that yields the Ernst potential from the functions a and 
/? can be solved for arbitrary coordinates (p, () ^ B, see pSM 7 ]. Moreover, 
the Ernst potential can be extended to the region B 7 since a, /3 G £> 7 . The 
different boundary values at the inner and outer sides of B follow from the 
construction. 

3. Finally, the proof of the third statement uses the linear system (|8]). For 
X G T with 3tX > 0, we establish the matrices $ ± [X, p = 0, ( = r B (l) ± 0] 
by integrating the linear system along the inner and outer side of B, with 
the initial values 



1 1 
1 -1 



The coordinate £ = r#(l) ±0 stands for the inner and outer side of B at 
p = 0, respectively. The integration of the linear system can be performed 
since for a sufficiently weak relativistic regular Ernst potential again the 
theorem by Picard and Lindeloff applies. At p = and ( = ± 0, 

we then establish $ from $ by $ = $M± where the regular matrices M± 
(defined on Y) are chosen such that $ assumes the structure @, i.e. we 
calculate the functions F± and B±. In order to determine 7 we choose some 
Sq > (say e = 1) and integrate the differential system flTT|) backwards, 
starting at e (here the initial conditions are just the established F±, B±), 
until we reach e — 0, and read 7 from the weak field expansion: 

7 = liml[(F_-F + ) + (£+-£_)]. 

So we have proven that the interior and exterior boundary data / and 
f,z j f,z (the derivatives also enter the linear system) of a sufficiently weak 
relativistic Ernst potential, which is regular within and without B, uniquely 
determine the function 7. But since the weakly relativistic Dirichlet bound- 
ary value problem is uniquely soluble [p4j| , the global interior and exterior 
Ernst potential (and hence the above derivatives) are determined by the 
boundary values / alone, which therefore determines 7 as well. 



7 Note that the linear integral equation might globally only be soluble for e < e* < e 
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Properties of the generalized hyperelliptic solutions. 

1. For ij; G M one finds 

/(e 2i ^£ 7 ) = /( £ 7)cos^ + isinV; 

cos ^ + i/(e7) sin ■?/> 

This can be seen by investigating /({X^jp+x , as P — ► oo with 



(12) 



H(X)i[(X-X v ), 



e 2i %(X) 



lim 7 P+ i = lim 

.R— >oo _R— >oo 



^(x-^)n(x-x, 

^=1 



lim 

ij— >oo 



^(X)(X-X p+1 )f](X-X, 



Note that in the integral terms of (||, [?D for z/ = p + 1 the substitution 
X = X p+ i(l + t 2 ) is useful. 

The property (|12|) describes the general invariance transformation of /(£ 7 ) 
which retains the asymptotic flatness (/ — > 1 as r — > oo). f(e 2l ^s r y) is 
obtained from /(ej) as one performs the transformation (21) of with 
the parameters a = sin -2 t/>, (3 = 7 = — cot r/> of that paper. 

2. A given function 7 : T — > C can be represented in many ways by sequences 
{{Xu and {if p }^. Consider the following example for a spherical 

boundary, = 1. Since 7 is analytic on T, it can be written as 

00 

7 (X) = (j^X*- 1 + lV X ~ 3 ) > X G r 



where 



7 



(±) 



2J-1 



Q± (cos $)Pj_i (cost?) sin^di? . 



Q±(cos$) = lim 



1 

2vri 



7 (X) rfX 
W z 



For each nGNwe find a representation 



7n (X) = ^( 7 f ] X^ 1 + TPX^') = P„e 2i ^X m " ]J(X - iW). 
3=1 



i/=i 



Thus we can write 

/(e 7 ) = lim f(ej n ) = lim 

with F n (X) = R n X m ". 



f({xl n) } Pn , eP ra ) cos^ n + i sin ip r , 
cosip n + if({xi n) } Pn , eH n ) sin ip n 
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3. If the analytic function 7 : T — > C is regular outside T and 7 — > as 
|X| — > 00, then one finds /(p, £; £7) = 1 within B and 



n . t -\\ \ G—s- sinh G 

/(p = 0,C>r B (l);67) = G cos h sinh G 

G cosh sinh G 

/(p = 0,C<-r B (-l);£ 7 ) = G- g , s inhG 



(13) 



with 



S± (X) = - ]^{X) ± 7 (X)J , G 2 (X) = s 2 + (X) - s 2 _(X), X without T. 

Hence, in the case of a flat interior field there is an explicit relation of the 
function 7 and the values of the Ernst potential along the exterior part of 
the symmetry axis. 

4. The Ernst potential / = f(e r y) possesses a reflectional symmetry 



t( n r.^A - ) f(P'C,£l) K for (p, C) within B , . 

J 1 /(P,C;^7)- K for (p, C) without 5 ' U4J 



if r B (r) = r B (— r) and 7(— X) = k7(J),k 2 = 1. We obtain the desired 
reflectional symmetry of the exterior Ernst potential for k — — 1. 

5. There is a weak field expansion of / = f(£j)' 

By differentiating (|| |^) with respect to e, where if and {X u } p are fixed, 
one deduces that 

i MrtW, , sE)] = ^ / H (X) f[(X - . 

The expansion (|i~5|) follows since, in the weak relativistic region, (|7|) yields 
X^ = X v + 0(e 2 ). 

6. Another interesting property is 



/( £T *) = /( £T ) for i*(X)= 7 (X). 

If, in particular, 7 = 7* then / = f{e^j) is real and belongs to the Weyl 
class. 
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3 Approximation of arbitrary exterior bound- 
ary value problems by Backlund type solu- 
tions 

Although the generalized hyperelliptic class permits the simultaneous solution of 
a sufficiently weak relativistic exterior and interior boundary value problem of 
the Ernst equation, its mathematical complexity makes it inconvenient for usage 
in a procedure to approximate the solution of the boundary value problem in 
question. Therefore one is led to investigate whether the much simpler Backlund 
type solutions suffice for our approximation scheme of only exterior boundary 
value problems. As we will demonstrate in this section, the freedom of the choice 
of the interior field allows us to set 7 7^ for X e T, which ensures that the 
corresponding Ernst potential / can be very well approximated by the Backlund 
type solutions. 



3.1 The class of the Backlund type solutions and its gen- 
eralization 

Definition. 

The Backlund type solutions G) depend on the set {^} ? of complex 

parameters and on the analytic function G : F — > C with G(X) = G(X) and are 
defined by 



where 



Mp,C; {Y u } q ,G) = f ^ 



D + = 



1 


1 


1 


1 


±1 


CK1A1 


a 2 A2 




1 


A? 


A2 " 


A 2q 


±1 


a 1 \l 


a 2 \\ 


a 2q^2q 



±1 ttxAf^ 1 a 2 \ 2 2 q - 1 ■ 
1 \\ q \ 2 2 q ■ 



X 2g-1 



' a 2q*2q 

,2q 
A 2q 



(16) 



fo 



1 r G(X) dX 



Y„-iz 



Ao;y A 



2v *2v-\ 
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. , A 2 „-i(r„ + iz) f G(X)dX \ 

= - tan H iS t w^xWJ ' 

The Ernst potential is only then regular inside and outside the 

boundary B if we additionally require that for each Y v G T there is G(YJ,) = 0.0 
For the desired procedure for approximating the solutions of exterior boundary 
value problems we will restrict ourselves to potentials fB{{Y u } q , G) with all Y v ^ 

r. 



Backlund type solutions as special hyperelliptic solutions. 



The Backlund type solutions form a subclass of the hyperelliptic solutions. In 
particular: 

f B ({Y v } q , G) = f({Y u Y u ...,Y q , Y q } , H) (17) 

with 

<? 

Y[{X-Y V )(X-Y V 



H(X) = G(X) 



and f({X u } p ,H) as defined in @. The proof of this works in the same manner 



as demonstrated in |23j for the solutions corresponding to disk like sources. The 
regularity of the solutions f({X u } p , H) at (p, () = (\5s[X v ]\, ^[X^]) applies in this 
specialization and is exhibited by the fact that oijXj is an even function in Xj. 
This means that / does not behave like a square root function near the critical 
points (p, C) — G^IX/] I; [^i/]) but rather like a rational function. 



Backlund transformations. 

The Ernst potentials (|T6|) are special Backlund transformations as described in 



53 (see formula 76) with complex conjugate parameters K 2u -i = Y„, K 2u = Y u 
and the real seed solution f = exp[U(G)]. The ctj in formula (76) of satisfy 
the Riccati equations 

\ /1 2\ fo,z 1 2\ fo,z 

a i,z = A^l - Qj ) — — , a jfg = — (1 - ocj) — . 

With constants Cj of integration, the general solution of these equations reads 

. fXjiKj + iz) I G{X)dX 

aj = - tanh n 3 , '- <p /v K ' + Cj 

3 \ 4vri J r {X - Kj)W z 3 

By the particular choice in fll6|) [i.e. C2V-1 = , C 2u = ivr/2], the atj become odd 
functions in Xj which ensures the regularity of the resulting Ernst potential at 

(p,C) = (l^]|^[n])- 

8 More precisely, if |G(y„)| > S > always holds, then the interior/exterior Ernst potential 
encounters a square root like behavior at (p,() = (I^P^IIj^P^]) when Y v tends to a point at 
the outer/inner side of T. 
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The generalization of the Backlund type solutions. 

The function 7 9 which belongs to the Ernst potential / = fB({Y u } q ,G) can be 
taken from flTTj): 

It is therefore a consequence of the above theorem that for given analytic functions 
G:r^CandS:r-^C with G(X) = G(X) and E(X) = -E(X), as well as 
sufficiently small e, the Ernst potential 

f B (~,sG)= \imf B ({Y^} q ,eG) 

q^oo 

exists and is independent of the particular choice of the sequence {{Y^^g}^ 
which serves to represent 5 by 

E(X) = lim In 

q^oo 

In particular we have 

f B (E,eG) = f(s 1 ) with 1 (X) = G(X)ex P [~(X)]. (18) 

From (|I8| ) we may conclude that any analytic function 7 which possesses a single 
zero X v e V and does not vanish at X v cannot be represented by two analytic 
functions G and E. Therefore, in the following we will restrict the discussion to 
7 ^ for X e r. 

Remark: As a particular consequence we now can confirm the conjectures 
formulated in ||23|| , sections 2 and 5. The first one of these conjectures treats the 
generalization of Backlund type solutions that describe disk-like sources of the 
gravitational field. In order to map the treatment of that place to our extended 
B here, we consider for given functions £ and g [] and an arbitrary point (p, Q 
without the disk, a boundary B which encompasses the disk but leaves the point 
(p, () outside B as well as all singularities of g(—X 2 /pl) and £(— X 2 / pfy outside 
T (p is the radius of the disk). If we now take any sequence {{Yv^^q}^ which 
serves to represent £ by 

£(x 2 ) = lim £,(* 2 ) , e 9 (x 2 ) = -ln 



9 In p3| . £ and g were assumed to be real analytic functions denned on the interval [0, 1]. g 
describes the seed solution while £ comprises the Backlund parameters. 



n 



X-Y 



(9) 



Y V (q) 



for x e Gq n 



irj g) -x 

■^iY^ + x 



n 



for x G [—1,1] 
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with all (pqYv^) being outside T, then the resulting Ernst potential /(/>, Cl £g> flO 
in at the chosen point (p, £) coincides with the above fs(p, C> — g> C<j) when 

G 9 (ip x) = 2(-l) 9+1 arccot(ix)#(x 2 ) , 
exp[S 9 (ip a;)] = (-l) q exp[-x£ q (x 2 )} . 

As g — > oo, the corresponding functions 7 g with 

Iqiipox) = G g (ip x) exp[E q (ip x)} 

converge towards the function 7, 

j(ipox) = — 2 arccot(ix) g(x 2 ) exp[— x^{x 2 )]. (19) 



So, due to the above theorem, the Ernst potentials /(p, Ci^flO m [HI tend to 
/(p, C)7)- Hence we can conclude, that the exterior gravitational field of a disk- 
like source is determined only by the local functions g and £ , and does not depend 
on a particular global representation of £ in terms of the parameters Y v . This is 
just the statement of the assumption in section 2. of |23|| . Note that the condition 
that Y v be outside the imaginary interval [— i, i] is always realized for sufficiently 
weak relativistic differentially rotating disks of dust since this is ensured by the 



positivity of the surface mass density a p (see appendix C.1.2 of [£3j; the function 
go is strictly positive and therefore the function £ is analytic). As one moves 
to more relativistic disks one might encounter a situation in which g and £ are 
not analytic (corresponding to the situation discussed above). However, for the 
examples investigated, this situation did not occur. 

In section 5. of ||23|| , conjectures were given regarding the functions g and £ for 
the disk-like hyperelliptic class of solutions. These conjectures are also proven 
by the above considerations. A specially chosen disk-like hyperelliptic solution 
as presented in appendix A of |^3| with an analytic function h and parameters 
{X u } p , coincides with /(p, (; 7) at a given point (p, () outside the disk, when 7 
reads 

p 

7(ip a;) = -2 arccot(ia;)7 Z )(ipoa;) , 7 D (ip a;) = h(x 2 ) \\(ix - X v ) 

u=l 

and a corresponding boundary B is taken in the above manner with the sin- 
gularities of 7d being outside Y. The comparison with ( pjf) yields exactly the 
conjectured functions g and £ in terms of h and {X u } p . It is of particular interest 
that thus the Neugebauer-Meinel-solution can be written as a well defined limit 
of Backlund type solutions. 

3.2 The approximation scheme and results 

In this section we describe our procedure to approximate the solution to an 
exterior boundary value problem of the Ernst equation. Due to the preceding 
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formulations, we seek an appropriate function 7 such that the corresponding 
Ernst potential / = fi^y) assumes prescribed exterior boundary values at B. 
As outlined above, we only consider functions 7 that do not vanish at T, and 
for which the corresponding Ernst potentials can be well approximated by the 
Backlund type solutions. We find in particular that in contrast to the Laplace 
equation, arbitrary boundary values of the Dirichlet type are not possible. 



The approximation scheme. 

1. At first we restrict the function 7 in an apropriate manner. Because of the 
desired reflectional symmetry (|14|) we take rs(r) = tb{— r) and 7(— A) = 
—7(A). To avoid zeros at T we demand ^(7) = const. = 70 7^ for A G T. 
This leads to 

7(A) = i 7o + A(A)+A 2 (A) 

with 



A 1 [X = r B (cos9)e 
A 2 [X = r B 



cosy e 



] = ^72fc-iCOs[(2A;-l)0] 

fc=i 

00 

] = ^72fcSin[2A;6>] , 7j e 



9 E [0,2tt]. 



fc=i 



The real constant 70 7^ may be chosen arbitrarily; the coefficients jj for 
j > 1 are then determined uniquely by the exterior boundary conditions. 

2. It is now possible to find analytic functions G and 5 with 

7(A) = i 7 (A)=iG(A)exp[S(A)], 

G(X) = 



A\{X) - [A 2 (X) + i 7o r , G [A = ir B (0)] = 70, 



H(A) 



arctan 



A 1 {X) + A 2 (X) 



7o 



+ arctan 



A 1 {X)-A 2 {X) 



7o 



+-ln 
4 



[A 1 (A) + A 2 (A)] 2 + 7q 2 
[A(A)-A 2 (A)f + 7 2 



[X = ir B (0)] = 0. 



With ([T2|), the Ernst potential reads then as follows 

f B &G)+i 



l + i/ fl (H,G) 



3. For a given analytic function H of the above kind we determine correspond- 
ing Backlund parameters {^} g in the same manner as outlined in ||23|| . 
This means that we solve the linear system: 



exp[H(^)]P(-^ 



-IYP(Z V 



u=l...q, 
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P(X) = f[(X-Y v ) = X* + fi i b J X* 

v=\ j=0 

to determine the coefficients bj and from these the zeros Y v of P. The q 
arbitrary different supporting points Z u G T are chosen to possess positive 
real and imaginary parts. 

4. By virtue of formula (16) we are now able to approximate the Ernst po- 



tential /(7) if the parameters {7^}^° are given. For the approximate nu- 
merical evaluation of {7j}? n from the given exterior boundary values we 
set 7j = for j > 2n and demand that the resulting exterior Ernst poten- 
tial 7(7) coincide with the prescribed boundary values at n different points 
{pj,(j) £ B with Cj > 0. This gives a complicated nonlinear set of In real 
equations to determine the unknown coefficients {7j}i n . As in the analo- 
gous treatment to solve boundary value problems for arbitrary differentially 
rotating disks of dust, we solve this system by means of a Newton- Raphson 
method. Again we provide good initial guesses by solving a sequence of 
boundary value problems with initially weak relativistic and finally the de- 
sired boundary values, and where the initial guess comes from the weak 
relativistic expansion (|Hi|). 

Results. 

The above approximation scheme has been executed for various Dirichlet bound- 
ary value problems. As expected from the above, there is always a weak rela- 
tivistic region within which the solution exists and can be well approximated by 
the Backlund type solutions. In particular, the exterior solution proves to be 
independent of the choice of the parameter 70, only the accuracy of the results is 
affected. 

As one moves further to more relativistic boundary values, there might be how- 
ever a limit beyond which the solutions cannot be extended. For example, take 
the exterior boundary values 



/[( P ,C)e £] = ! + £ 



3 ,2 
— + C 
2 s 



with the spherical boundary B characterized by = 1. One finds a limiting 
parameter e ~ 0.68 beyond which the boundary value problem does not seem 
to be soluble. This has been verified not only by various choices of the above 
parameter 70 but also by a two-dimensional numerical method, which solves the 
exterior Ernst equation directly without using analytic solutions^. All these 
routines yield the same limiting parameter Eq which is a strong indication that 



10 For this method, we insert a two-dimensional Chebyshev expansion of the exterior Ernst 
potential into the Ernst equation and find the Chebyshev coefficients. 
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for e > £q the Ernst potentials corresponding to the above boundary values do 
not exist. 
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